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ABASTRACT 

A  micromechanical  model  for  three-dimensional  open-cel!  foams  is  developed  by  using  the  matrix 
method  for  space  frames  in  structural  mechanics  and  tetrakai decahedral  unit  cells,  The  effective  elastic 
properties  of  foams  are  determined  employing  unit  cells  subjected  to  three  different  inodes  of  loading. 
The  thirty-six  struts  of  each  tetrakaidecahodral  unit  cell  are  treated  as  uniform  slender  beams,  and  the 
twenty-four  vertices  as  rigid  joints.  All  four  deformation  mechanisms  (i.e..  stretching,  shearing,  bending 
and  twisting)  possible  under  the  specified  loadings  are  incorporated,  and  four  different  stmt  cross  section 
shapes  (i.e,,  circle,  square,  equilateral  triangle  and  Plateau  border)  are  treated  in  a  unified  manner.  The 
formulas  for  determining  the  effective  Young’s  modulus,  Poisson's  ratio  and  shear  modulus  of  open-cell 
foams  that  are  undergoing  linearly  elastic  deformations  are  derived  using  the  composite  homogenization 
theory,  The  new  formulas  include  all  necessary  parameters,  unlike  those  provided  by  existing  models. 
These  formulas  indicate  that  the  foam  elastic  properties  depend  on  the  relative  foam  density,  the  shape 
and  size  of  the  strut  cross  section,  and  the  Young’s  modulus  and  Poisson’s  ratio  of  the  strut  material.  By 
applying  the  new  model,  a  parametric  study  of  sample  cases  is  conducted  for  carbon  foams,  whose 
modeling  motivated  this  study.  The  predicted  values  of  the  effective  properties  agree  favorably  with  those 
based  on  existing  models  and  experimental  data  for  the  Mode  1  loading  case,  which  is  the  only  case  that 
has  been  well  studied  in  the  literature.  Comparisons  of  the  effective  clastic  properties  for  the  three  loading 
cases  quantitatively  show  that  carbon  foams  exhibit  certain  degree  of  anisotropy. 

KEYWORDS:  Micromechanics;  Open-celi  foams;  Effective  elastic  properties;  Matrix  method; 
Composite  homogenization  theory 

1.  INTRODUCTION 

Microcellular  graphitic  carbon  foams  first  developed  at  the  Air  Force  Research  Laboratory  (AFRL)  in 
1990s  (e  g..  Hall  and  Hager,  1996)  are  rapidly  emerging  as  a  new  class  of  ultra-light  cellular  materials  for 
structural  and  thermal  management  applications  because  of  their  excellent  mechanical  and  thermal 
properties  (Roy  et  al,,  1998).  Reliable  structural  applications  of  these  carbon  foams  depend  on  accurate 
understanding  of  their  mechanical  behavior.  Carbon  foams,  like  other  cellular  solids,  are  topology- 
sensitive  and,  therefore,  their  mechanical  modeling  requires  the  incorporation  of  cell’s  micrastructura], 
features, 

AFRL  graphitic  carbon  foams  are  blown  from  anisotropic  pitch  through  a  bubble  forming  process  (e.g,, 
Anderson  et  al.,  2000),  and,  as  a  result,  microstructures  of  the  solidified  carbon  foams  are  controlled  by 
the  principle  of  minimum  surface  energy.  As  shown  in  Fig.  ! ,  the  mic restructure  of  an  AFRL  carbon 
foam  has  a  three-dimensional  (3-D)  open-ceil  topology,  As  a  first  approximation,  a  regular 
tetrakai  decahedron  (see  Fig.  2(a))  can  be  adopted  as  the  repeating  unit  to  represent  this  foam 
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microstructure,  Tetrakai  decahedron  i$  known  to  be  (he  only  polyhedron  that  can  pack  with  identical  units 
to  Till  space  and  nearly  minimize  the  surface  energy  (e,g.f  Weaire  and  Fortes,  1994).  With  all  of  its 
vertices  connected  by  struts  and  each  vertex  shared  by  four  struts,  such  a  fourteen -sided  polyhedron  can 
be  generated  by  uniformly  truncating  the  six  comers  of  an  octahedron  and  contains  eight  regular 
hexagonal  faces  and  six  square  faces. 


Fig,  1 .  Micrograph  of  an  AFRL  carbon  foam  (Roy  et  til.,  1 998), 


(a)  (b)  (c) 


Fig.  2,  Tetrakaidecahedral  unit  cell  under  loading:  (a)  Mode  I,  (b)  Mode  II,  (c)  Mode  ID, 

Several  micromechanical  models  for  3-D  open-cell  foams  have  been  developed  utilizing 
tetrakaidecahedral  unit  cells,  which  are  reviewed  in  Li  et  <il.  (2003).  Very  recently,  a  mteromechanics 
model  for  predicting  effective  Young's  modulus  and  Poisson's  ratio  of  such  foams  was  provided  by  Li  et 
al.  (2003)  using  an  energy  method  based  on  Castigliano’s  second  theorem,  whose  predictions  compare 
favorably  with  those  based  on  available  experimental  data  and  Finite  element  analyses.  However,  only 
Mode  1  loading  (see  Fig,  2)  was  considered  there,  as  was  the  case  with  almost  all  previous  studies  on  the 
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topic,  Also,  how  to  predict  effective  shear  modulus  was  not  addressed  in  Li  et  al.  (200,1).  When  loads  are 
not  applied  on  two  opposite  square  faces  of  a  tetraka [decahedral  cell,  the  elastic  properties  obtained  might 
be  different  from  those  when  the  cell  is  subjected  to  Mode  [  loading.  For  example,  cubic  symmetry  is 
present  under  Mode  I  loading,  but  is  absent  under  Mode  II  or  III  loading.  For  the  latter  cases,  anisotropy 
involved  can  be  significant.  Hence,  new  micromechanicai  models  capable  of  accounting  for  the 
differences  in  mechanical  behavior  of  foams  (cells)  under  three  distinct  modes  ofloading  are  necessitated. 
The  objective  of  this  paper  is  to  present  such  a  new  model.  The  tetrakaidecahedral  unit  cells  illustrated 
in  Fig.  2  are  adopted,  and  the  matrix  method  for  space  frame  structures  and  the  homogenization  theory  for 
composite  materials  are  employed  in  the  formulation. 

2.  MATRIX  METHOD  FOR  SPACE  FRAME  STRUCTURES 


2.1  Method 

In  this  matrix  method  the  members  of  a  space  frame  structure  are  regarded  as  rigidly  connected  at 
joints  (or  nodes),  each  of  which  has  six  degrees  of  freedom  -  three  translational  and  three  rotational  - 
when  unrestrained.  The  twelve  possible  displacements  of  the  two  joints  of  a  member  may  be  described 
relative  to  the  (global)  reference  axes  x.  y  and  z  or  to  the  (local)  member-oriented  axes  yM  and  Zu,  as 
shown  in  Fig.  3.  The  xu  axis  is  set  to  be  along  the  axis  of  the  member  from  i  to  j,  and  the  yw  and  zu  axes 
are  chosen  as  the  principal  axes  of  the  cross  section  at  joint  i  of  the  member, 

In  terms  of  the  member  coordinates  [jt.Wi  y.M,  z*t,  I,  the  forces  at  the  ends  of  a  member  are  related  to  the 
displacements  at  the  ends  by  the  member  stiffness  matrix,  i.e,, 


where  and  Sji  are,  respectively,  displacement  vectors  at  nodes  i  andy,  and  and  fy,  respectively,  force 
vectors  at  nodes  i  andy  for  member  f\  The  member  stiffness  matrix  in  Eq,  (1)  is  composed  of  four 
sub  matrices  kf,,  k^,  and  k^ ,  The  first  subscript  of  each  submatrix  represents  the  forces  arising  at 

the  end  associated  with  this  subscript  caused  by  unit  displacements  occurring  at  the  other  end  (associated 
with  the  second  subscript). 


When  expressed  in  terms  of  the  reference  coordinates  {jr,  y,  ;|,  Eq.  (I)  becomes 


(2) 
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where  Fy  and  Fji  are,  respectively,  force  vectors  at  nodes  i  and  A|  and  Aj  are,  respectively,  displacement 
vectors  at  nodes  i  and  j,  and 


K«  ”T„-,ki'T,,Ku  =Vk„T#,  KJ(  .T,-'k,T,.  =T;‘k^T„, 


(3a-d) 


are  the  global  stiffness  submatrices,  with  being  the  transformation  matrix  for  member  “i,  f  (from 
global  to  local). 

Nodal  equilibrium  requires  that  the  applied  loads  at  a  node  be  balanced  by  the  internal  member  end 

forces.  Suppose  that  members  “t,  a'\  '*t,  b" . "j.  n "  are  connected  to  node  /.  Then,  the  equilibrium  of 

node  i  requires 


P)  ”P|»  +Flb  +  ’  +  Pin’ 

where  P|  is  the  applied  load  vector  at  node  i.  Using  6q,  (2)  in  Eq,  (4)  gives 

P|  =  KiA  +  K|bAh  +  ",  +  Kl0A11, 

where  K„  =  K  J  +  K  J  +  *  ■  ■  +  K  J .  Subsequently,  enforcing  the  equilibrium  of  al  i  nodes  results  in 


(4) 

(5) 


P, 


K.i  K, 

KJ(  K 


K 


IN 


a 


‘IN 


Nl 


K 


N2 


K 


NN 


A, 

A, 


(6) 


In  Eq.  (6),  which  is  called  the  initial  stiffness  equation,  each  nodal  force/displacement  vector  has  six 
components,  The  coefficient  matrix  of  these  equations,  denoted  as  Ki,  is  known  as  the  initial  structure 
stiffness  matrix.  Each  submatrix  of  K|  is  a  6  x  6  square  matrix.  Equation  (6)  can  be  rearranged  as 


Ifr  1  ^  fp  j  K  FR  |  Af  | 

iM’LkViKoKr 


(7) 


where  the  subscripts  F  and  R  refer  to  free  and  restrained  displacements,  respectively.  From  Eq,  (7)  the 
free  nodal  displacements  can  then  be  obtained  as 
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2,2  Member  Stiffness  Matrix 

By  following  the  procedure  used  in  Spillers  (1972)  for  2-D  frames,  the  member  stiffness  submatrices 
for  3-D  frames  undergoing  transverse  shearing  deformations  as  well  as  axial,  flexural  and  torsional 
deformations  are  obtained,  for  the  first  time,  as 
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where  d,  L  are,  respectively,  the  cross-sectional  area  and  length  of  the  member,  /T  and  /,  are  the  second 
moments  of  area  of  the  member  cross  section,  J  is  the  polar  second  moment  of  area,  £  and  G  are, 
respectively,  the  Young’s  modulus  and  shear  modulus  of  the  member  material,  and  kf  and  fc,  are  the 
transverse  shear  factors.  It  is  noted  that  when  ft,  and  it,  approach  infinity,  Eqs.  (9a-d)  will  be  reduced  to 
the  member  stiffness  sub  matrices  for  members  involving  only  axial,  flexural  and  torsional  deformations 
{e,g.,  Balfour,  1986),  Substituting  Eqs,  (9a-d)  for  the  member  stiffness  submatrices  into  Eqs.  (3a-^d) 
results  in  the  following  global  stiffness  submatrices: 


(Kjffi=T-'(kJ^,  (K^T-’tk^TV  <I«a-d> 


3.  FORMULATION 

Analysis  of  a  space  frame  structure  can  be  conducted  by  following  a  computer-oriented  procedure 
outlined  here.  First,  the  6  x  6  global  stiffness  submatrices  in  terms  of  the  reference  coordinates  |jc,  y,  z) 
are  generated  for  each  member  in  the  structure  (e.g„  member  ‘7,  /’  shown  in  Fig.  3)  using  Eqs.  ( lOa-d). 
Member  “i,  j”  contributes  to  the  stiffness  matrices  of  nodes  i  and  j  at  the  ends  of  the  member,  Hence, 
corresponding  elements  from  the  global  stiffness  submatrices  for  this  member  may  be  transferred  to  the 
overall  stiffness  matrix  by  use  of  suitabie  indices  (e.g,,  Weaver  and  Gere,  1990),  The  six  possible 
displacements  at  a  particular  node  i  are  denoted  by  the  following  indices: 

6(i  -  1)  +  I  =  index  for  the  x  component  of  translation, 

6(i  -  1 )  +  2  -  index  for  the  y  component  of  translation, 

6(f  -  ! )  +  3  =  index  for  the  z  component  of  translation, 
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6(#  -  1)  +  4  =  index  for  the  x  component  of  rotation, 

6(1  -  1 )  +  5  =  index  for  the  v  component  of  rotation, 

6{i  -  S)  +  6  =  index  for  the  z  component  of  rotation. 

The  transferring  ni!e  is  that  element  (i'iJ,)  (i,  J,  «  {1,2,  ■ 6J)  of  submatrices  (k*  }. .  (kJ^,  (kJ 

and  {k^  is  added,  respectively,  to  locations  (6i  -6  +  iu  &  -  6  +  j,),  (6/  -  6  +  iu  6J-6+  jx),  (6/  -  6  +  ilt 

6i  -  6  +  ji),  (6j  -  6  +  iu  6/  -  6  +  jt)  in  the  initial  structure  stiffness  matrix  Kj.  Construction  of  the 
complete  Ki  involves  forming  and  transferring  (k^,  (kJ^  .  (Kj(  )m  and  (KjJ/or  all  members  of  the 

structure. 

The  second  phase  of  the  analysis  is  to  generate  initial  load  vectors.  As  shown  in  Fig,  3,  the  actions 
applied  at  node  /„  (P|)i,  (PJj  and  (Pj)i,  are  the  x,  y  and  z  components  of  the  concentrated  forces,  and  the 
actions  (Pi)*  (PJs  and  (PJ4  are  the  x,  y  and  ;  components  of  the  moments  applied  at  node  f.  They  are  put, 
respectively,  at  the  locations  (Si  -  5),  (6i  -  4).  (6 i  -  3),  (Si  -  2).  (6i  -  I)  and  (60  in  the  initial  load  vector. 

After  the  construction  of  the  required  matrices  is  completed,  substituting  them  into  Eqs.  (7)  and  (8) 
finally  gives,  with  the  boundary  conditions  enforced,  the  solution  for  free  nodal  displacements. 


3.1  Mode  I  Loading 

The  tetraka [decahedral  unit  cell  under  Mode  I  loading  is  shown  in  Fig.  4,  with  all  of  its  members  and 
joints  numbered.  All  of  the  members  (struts)  of  the  cell  have  the  same  length  L  and  cross-sectional  area  A. 
In  addition,  the  following  relations  can  be  identified  (Li  and  Gao,  2003)', 

2  Jo 

A  =  —  RL\  /r  =  /(=cAJ,  J  =c,AI,  (11) 

where  f?  is  the  relative  foam  density,  and  c  and  C|  are  two  geometric  constants  depending  on  the  strut 
cross  section  shape. 

By  following  the  procedure  outlined  above  and  using  the  structural  and  loading  data  illustrated  in  Fig. 
4,  the  nodal  displacements  for  each  non-restrained  joint  of  the  cellular  structure  can  be  obtained.  The 
effective  Young's  modulus,  Er,  of  the  unit  cell,  and  thus  of  the  Foam,  can  then  be  determined,  using  the 
average  strain  theorem  In  the  homogenization  theory  (e,g„  Gao  et  al.,  2003),  as  (Li  and  Gao,  2003) 


£'  = 


2j2 P 
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(12) 


where  A,  [=  (A  1)2]  is  the  displacement  of  the  unit  cell  in  the  negative-y  direction.  From  Eqs.  (9a-d)  and 
(1 1)  it  follows  that  that  A*,  and  thus  Ef ,  depends  on  the  relative  foam  density  (R),  the  size  and  shape  of 
the  strut  cross  section  (c.  cx,  ky  and  fc.),  and  the  elastic  properties  of  the  strut  material  (£  and  G), 

The  effective  Poisson’s  ratio  u^.can  be  obtained  as  (Gao  and  Li,  2002;  Li  and  Gao,  2003) 


*  8  A,  y  40 


(13) 


where  (At)i  stands  for  the  ^-component  of  the  (translational)  displace  men  t  a  I  node  i,  and  jr(0  denotes  the 
x-eoordinate  of  node  i.  Clearly,  1/  =  v'v  due  to  symmetry. 


Finally,  the  effective  shear  modulus  G^can  be  shown  to  be  (Li  and  Gao,  2003) 

9V2  P 
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where  A„  and  Aw  are,  respectively,  the  displacements  of  the  cell  in  the  xs-  and  ^-directions,  which  are 
oriented  45°  relative  to  x-  and  y-axes  (see  Fig,  5). 


Fig,  4,  Tetrakaidecahedral  unit  cell  under  Mode  I  loading. 


Fig.  5.  Equivalent  loading  configuration  for 
shear  modulus  (Mode  I  loading). 


Fig.  6.  Configuration  for  Young's  modulus 
and  Poisson's  ratios  (Mode  fl  loading), 


3-2  Mode  11  and  Mode  Hi  Loading 

Formulas  for  predicting  the  effective  foam  properties  in  the  Modes  II  and  III  loading  cases  can 
be  derived  by  following  procedures  similar  to  that  used  above  for  the  Mode  I  loading  case. 

For  the  Mode  II  loading  case  (Fig.  6),  the  effective  Young’s  modulus  is  obtained  as  (Li  and 
Gao,  2003) 


9^2  P 
4LA’t  ' 
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and  the  effective  Poisson's  ratios  as 


■  _  ^  (A, ),  *  _  V*6  L  (At)5 

20  A'l  f  x(i)  ’  iT“  24  A"r^C0" 


(lfia,b) 


where  A”  and  A“  are,  respectively,  the  displacements  in  the  loading  and  lateral  directions  of  the 

unit  cell  under  Mode  [Noading,  z{i)  denotes  the  ;-coordinate  of  node  t,  and  (A,))  stands  for  the  z- 
component  of  the  {translational)  displacement  ai  node  L 

For  the  Mode  111  loading  case  {Fig.  2c).  the  effective  Young's  modulus  is  obtained  as  (Li  and 
Gao,  2003 ) 


5^2  P 

£'=iZS 


Hi 


and  the  effective  Poisson’s  ratios  as 


22  A1"  ,  jc{f)  ’  ^  16  A™  i  z(i) 


(17) 


(I8a,b) 


where  A”1  and  A1’1  are>  respectively,  the  displacements  in  the  loading  and  lateral  directions  of  Lhe  unit 
ce!]  under  Mode  III  loading. 


4,  NUMERICAL  RESULTS 


To  illustrate  the  new  model  a  parametric  study  of  sample  cases  is  conducted  here  for  three-dimensional 
open-cell  carbon  foams,  whose  modeling  motivated  the  present  study.  The  M ATLAS  program  (of  The 
Mathworks,  Inc.)  is  used  in  the  computations.  Young's  modulus  (£)  and  Poisson's  ratio  ( v)  of  the  carbon 
strut  material  which  is  assumed  as  isotropic,  are,  respectively,  taken  to  be  15.61  GPa  and  0-33.  and  the 
maximum  value  of  the  relative  foam  density  (R)  to  be  0.22,  as  was  done  in  Li  et  ai.  (2003).  The  shapes 
used  here  for  the  strut  cross  section  also  include  circle,  square,  equilateral  triangle  and  Plateau  border. 

The  predicted  effective  Young’s  modulus  (Ey\  Poisson's  ratio  (v„*)  and  shear  modulus  (Gzy )  varying 

with  the  relative  density  (R)  for  carbon  foams  with  Plateau  border  strut  cross  sections  and  under  Mode  I 
loading  are  illustrated  in  Figs.  7+  S  and  9,  where  they  are  also  compared  to  she  predictions  of  several 
existing  models. 
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Fig.  7.  Young's  modulus  vs.  re  I  alive  density.  Fig  8.  Poisson's  ratio  vs.  relative  density. 

From  Fig,  7  it  is  seen  that  the  values  of  Ey  predicted  by  the  current  model  and  those  by  the  four  earlier 
models  are  very  close  for  low-density  foams  (with  R  <  0.1).  Figure  8  shows  that  decreases 
monoton  tea  by  as  R  increases  according  to  the  new  model  This  trend  is  the  same  as  that  predicted  by 
three  other  models.  Also,  values  of  predicted  by  the  current  model  (e.g.,  v^*  =  0.346  at  R  =  0,1)  arc 
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dose  to  1/3,  which  is  a  value  suggested  by  Gibson  and  Ashby  (1997)  based  on  available  experimental 
data  (ra  n  gi  ng  from  0. 1 5  to  0 . 4 ). 
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Fig.  9.  Shear  modulus  vs.  relative  density. 


Fig.  10-  Young's  modulus  vs.  relative  density. 


It  is  dear  from  Fig.  9  that  the  values  of  G*y  predicted  by  the  current  model  are  closer  to  those  provided 
by  Gibson  and  Ashby  (1997)  based  on  the  experimental  data  than  the  values  given  in  Warren  and  Kraynik 
(1988).  Moreover,  Fig.  9  reveals  that  rhe  predictions  of  GZf  by  the  three  models  are  very  dose  for  low- 
density  foams  (with  R  <  0.1), 

Figure  10  compares  the  predicted  values  of  the  effective  Young’s  modulus  (  Ey )  of  carbon  foams  with 
Plateau  border  strut  cross  sections  when  the  unit  cell  is  under  Modes  I,  II  or  HI  loading.  It  is  shown  that 
Eif  increases  monotonically  for  all  the  three  loading  cases  as  R  rises.  The  values  of  for  the  Mode  l 

and  H  loading  cases  are  very  close  for  low-density  foams  (with  R  <  0. 1 5)t  which  conforms  to  the  weak- 
anisotropy  observed  by  Warren  and  Kraynik  (1997).  In  this  same  low-density  region,  however,  the  values 

of  E*  for  the  mode  ID  loading  case  are  50%  lower  than  those  for  the  other  two  loading  cases.  This 

reflects  the  degree  of  anisotropy  of  carbon  foams  as  measured  by  £y * 

For  the  tetrakaidecahedra!  unit  cell  under  Modes  I,  II  or  III  loading,  the  predicted  values  of  the 
effective  Poisson's  ratios  varying  with  R  are  illustrated  in  Figs.  1 1  and  12.  For  the  Mode  II  loading  case, 

v\^  is  0.627  at  R  =  0  and  decreases  gradually  with  the  increase  of  Rr  while  ^rv  starts  to  descend  from 
0,212,  When  the  unit  cell  is  loaded  at  two  farthest  joints  (Mode  III  loading),  the  values  of  both  and 
v[v  are  very  close  to  those  for  the  Mode  II  loading  case.  However,  there  are  large  differences  between 

vxr  and  v\.  for  the  same  loading  case.  Furthermore,  Figs.  II  and  12  show  that  for  V ^  or  the 

differences  between  the  values  for  the  Mode  I  loading  case  and  those  for  the  Modes  E  and  III  loading 
cases  are  significant.  This  indicates  that  carbon  foams  do  exhibit  distinct  anisotropy  when  viewed  in  terms 

of  Poisson's  ratios-  Also,  it  is  observed  from  Fig,  12  that  the  predicted  values  of  v]y  for  all  of  the  three 

loading  cases  fall  within  the  range  of  0.15  -  0.4,  which  is  the  experimental  data  range  quoted  in  Gibson 

and  Ashby  (1997). 

The  predicted  effects  of  the  relative  foam  density  [R]  on  the  effective  shear  moduli  ( G ^  and  Glf)  are 

illustrated  in  Figs.  13  and  14.  It  is  dear  that  the  values  of  shear  moduli  are  very  close  for  all  the  three 
different  loading  cases.  This  agrees  with  the  similarity  between  the  results  of  the  effective  shear  moduli 
for  the  Modes  I  and  II  loading  cases  reported  in  Warren  and  Kraynik  (1997).  In  addition,  it  is  noticed  that 

the  differences  between  the  values  of  Gr.  and  G\  under  the  same  mode  of  loading  are  indistinguishablc- 

These  imply  that  the  anisotropy  as  measured  by  the  shear  moduli  is  insignificant. 
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Fig.  11,  Poisson's  raiio  vs.  relative  density. 
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Fig.  13.  Shear  modulus  vs.  relative  density. 


Fig.  14.  Shear  modulus  vs,  relative  density. 


Finally,  the  order  of  influence  of  strut  cross  section  shapes  on  the  effective  properties  of  carbon  foams 
appears  to  be  the  same  in  all  three  loading  cases.  That  is.  Plateau  border,  equilateral  triangle,  square  and 

circle  strut  cross  sections  form  a  descending  order  for  £y,  Giy  and  G,v ,  and  a  ascending  order  for  both 

vx,  and  vty .  At  R  =  0, 1,  the  maximum  relative  differences  of  the  three  clastic  moduli  and  two  Poisson's 

ratios,  resulting  from  the  use  of  different  shapes  of  strut  cross  section,  for  the  three  loading  cases  are 
summarized  in  Table  1 ,  In  all  of  the  three  loading  cases,  shear  moduli  are  most  sensitive  to  the  strut  cross 
section  shape,  which  is  followed  by  Young's  modulus  and  Poisson’s  ratios.  It  is  noted  from  Table  1  that 
the  values  of  relative  differences  between  different  loading  cases  are  very  close  for  all  five  elastic 
properties  except  tor  Poisson's  ratios  in  the  Mode  III  loading  case,  which  are  substantially  lower  than 
those  in  the  Modes  I  and  II  loading  cases.  This  supplements  the  observations  about  anisotropy  in  foam 
elastic  behavior  drawn  from  Figs.  1 1  and  12. 


Table  I  Maximum  relative  differences  due  to  the  use  of  different  strut  cross  section  shapes 
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